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ABSTRACT 


The purpose of this thesis is to provide the student 
of queuing systems with a vehicle through which a better 
understanding of the interrelationships between stochastic 


processes and queuing systems can be achieved. 
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is ENERODUGTION 


A queuing system is comprised of action facilities 
(servers) and waiting facilities (queues or lines). The 
relationship between the components can be of many forms, 
esc with their own characteristics. Consider a tool 
€rib in’ a machine shop where the machinists are obliged 
to check out their Eagles or iconsider the itelephone system, 
in which callers request communication capacity. Both 
Situations involve requests for action from limited resources, 
iecnese cases, servers or telephone circuits. . In order to 
a@ibacate Gesources properly, it is necessary to understand 
the interaction between requests for action and the actions 
themselves. Individually these can be modeled as stochastic 
processes, but in order to predict flow through the queuing 
Syorem, 1t 1S necessary to understand the interrelation 
between the patterns of event occurrences. When demands 
occur at regular intervals and actions are of fixed dura- 
Puon, it is relatively easy to predict flow; if there are 
more requests per unit of time than actions over the same 
time span, the queue will grow without bound unless the 
event relationship is changed. When variation is added to 
the timing of events, the flow is not as obvious. For 
example, if the average number of requests for action 
equals the average number of actions, per unit of time, but 
both Occur Wuth variability, then it would be useful to be 


able to predict such quantities as expected delay and action 


— + oe 


facilaty sdle time. Queuing models are used for’ this 
purpose. 

The purpose of this thesis is to provide the student 
Of queuing systems with a vehicle through which a better 
understanding of the interrelationships between stochastic 
processes and queuing systems can be achieved. 

This is done through a combination of a computer simula- 


tion and a computer-driven graphical display. 


II. QUEUING MODEL SIMULATION 


This computer program simulates queuing systems in two 
Segees.. First, a next-event type computer simulation creates 
a time history of the queuing system. Second, a dynamic 
Sraphical display depicts the occurrence of the queuing 
process in real time. 

Five basic, academically interesting, queuing models, 
each with parameter options, can be simulated. Queuing 
models are generally described in terms of the distribution 
Gi request occurrences (arrival distribution), the distribu- 
t2onm of action occurrences (service distribution), and the 
number of action facilities (servers). Specifically, the 
times between request occurrences are taken to be indepen- 
dent and identically distributed realizations from the 
aenaval distribution. The duration of each action occur- 
Hence 1S also taken to be an independent and identically 
distributed realization from the service distribution. 
Variations on these assumptions occur for some models and 
are described in Section A. 

The abbreviation G/G/C describes a queuing system in 
Which the arrival and service distributions are general 
(unspecified) and there are C servers. 

ihe! five options in. this simulation are: 

G2 G/G/C 
(2) G/G/C with losses 
(3) G/G/C with feedback 


(4) G/G/C with finite source 

(S) Two simultaneous G/G/C queues for comparison. 
The specific choices of arrival distribution and service 
distribution which can be simulated are fully described 


im’ Sect ien, Be 


A. THE QUEUING MODELS 
1. G/G/E 


Prcturer 1 


This is the simplest of the queuing models considered 
here. Each arrival waits in the queue until a server is 
available, enters service for a calculated period of time, 
then departs the system. 

This model may be simulated with from one to six 


SETVEES:. 
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2. G/G/C.with, Losses 


Picture 2 


Consider making a long distance telephone call on 
Mother's Day. If all the trunk lines are busy, your request 
for service will be rejected. The G/G/C model with losses 
hemlects this by rejecting arrivals who find the queue 
Pei. «The queue. limit used is, specified as an input 
parameter. Note that rejected arrivals are lost and are 
met assumed to return later. The fraction of all arrivals 
who are served (not rejected) is recorded for subsequent 


display. 
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Picture 3 


Consider a closed-circuit educational television 
myouem ate da School. The finite population consists of all 
the televisions in the system. When one fails, it enters 
fie queue to be Serviced. After repair, it re-enters the 
population pending the next failure. 

This situation is modeled by assuming that each 
action source, upon entering the idle population, is 
assigned a waiting time from the arrival distribution. 
When this time expares, the source requests service, 
joanine the queue 1f mecessary. As each action source 
terminates its service, it returns to the idle population. 
The user specifies the number of servers (one to six) and 


che "size of the population. 
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4. G/G/C with Feedback to the Queue 


Picture 4 


This is an extension of the G/G/C model, with from 
one to six servers, allowing customers leaving service to 
depart the system or to return to’ the’ queue. ‘At each 
service completion, the customer leaves the system with 
probability P’ or’ returns* to the* queue’ with probability 1-P. 


The user specifies P as an input option. 
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9. Two Simultaneous G/G/C Queues 


Picture 5 


This option allows comparison of two standard 
G/G/C queues, each with up to three servers. All parameters 
of each system are specified separately, and, if desired, 
identical arrival streams may be utilized. This allows 
comparisons between different arrival and/or service distri- 
butions and associated parameters, and between a different 
number of servers. Also, the effect of the various queuing 


disciplines can be compared. 


B. ARRIVAL/SERVICE DISTRIBUTIONS AND QUEUE DISCIPLINES 
As well as options already described, the following 
arrival and service distributions are available to the 


user of the program. Any of the following probability 
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aTstrabutLons Can be uSed as -an arrival or service 
a@rstribution. 
i DiSteriDuUcLons 


a. K - ERLANG (exponential when K=1) 


Ke J 
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Input parameters are: 
(1) kK 


(2) Rate - in terms of mean number of events 
PEL Manwce:. 


b. Hyperexponential 


BCC) = Pelee 1c) - ae 2 


where: 


Ga pee qe =*4h 


(2) «eV = 
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Input,;parameters ane: 

(1) Rate, specified in terms of mean number 
of events per minute. 

(2) Coefficient of. variation,» (must. be- greater 
than, 220) 


e¢. . Degenerate (Fixed: Rate) 


Ges Ne 


Input parameter: 
(1) Rate A specified in terms of mean number 
of events per minute. 
2. Queuing Disciplines 
For comparison purposes, there are four queuing 
disciplines to choose from. These determine the priority 
Byewhich members of the queue will enter service. 
| as -FIFO = the £irst* anto queue: is the’ first: out. 
Over all disciplines which are independent of service. times, 
this one minimizes the variance of delay. 
b. LIFO - the last into the queue is the first 
out. Some systems, by their nature, function in this manner. 
c. SSTF - the queue member with shortest service 


time is served first; this minimizes the mean delay over 
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wi other disciplines, but the vance of delay is greater 
taan that of FIFO. 

d. RANDOM DRAW - applies to a system in which 
insufficient information is known to implement any other 
discipline. 

3. Length of the Simulation 

The period of simulation is specified in minutes 

Byethe user. All queues are empty and servers idle at the 


seart Of the simulation. 
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Iii. THE METHOD 


A. THE COMPUTER SIMULATION 

The computer simulation is a modified next event type. 
There are two types of events, arrivals to the queue and a 
completion of service, With the first there will be an 
associated entry to service if’ there 15 an idle server, 
and with the last an associated entry to service if the 
queue is-not empty. 

The simulation is begun with the generation of the first 
amiaval time and 1ts- associated service time. The times are 
generated via the transformation of random numbers to stochas- 
tee variates of the type specified. At the occurrence of 
an arrival, the next arrival time and associated service 
time are generated. When a server is available, a waiting 
member of the queue enters service with completion time 
equal to the sum of his assigned service time and the time 
hemdeparts the queue. Thus with the clock starting at zero, 
a pine history of events is created which tells what happened 
and when. 

At the same time, the appropriate values are retained 
from. which statistics of the simulation can later be 


calculated. 


B. THE GRAPHICAL DISPLAY 
At the completion of the simulation phase, the time 


history is depicted dynamically on the cathode ray tube 
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timenedl time. Pictures: t to 5 show the form of the display , 
with the man-like images in motion along the appropriate 
path, and the number in the queue box showing the queue's 
current length. 

The timing and motion of the display are controlled 
by the analog computer via the hybrid connections to the 
digital computer. The first connection is for the purpose 
Ge timing Until the next, event iS to be =initaated.- In this 
case, the analog computer integrates from zero to a preset 
came value in volts, at the rate of one volt per second. 
ihe interrupt to the appropriate subroutine then occurs, 
and the analog computer is reset with the next time value, 
after which the integration is reinitiated. 

The second connection causes the dynamic motion to 
eecur by continually checking events’ in progress and causing 
the position of the corresponding images to be incremented 
appropriately; this interrogation and incrementation occurs 
approximately 10 times: per second depending on the settings 
6 the analog computer. Movement increments are between 
me -halt inch and two inches (depending on the path), 
causing a motion of between five inches and 20 inches per 


second. 
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IV. STATISTICS FROM THE SIMULATION 


STATISTICS 
02 


3.99 
.95 
94 


a 
.63 
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ENTER 1. CONTINUE ‘ 
2. PARAMETER LISTING 


At the end of the simulation period, the following 
statistics are calculated and displayed, as shown above. 
AT = time - average length of the queue 
BT = time - average number of busy servers 


LT = time - average number in the system = QT + BT 


QC = customer - average length of the queue 
BC = customer - average number of busy servers 


LG = customer - average number in the system = QC + BC 


DE = customer - average delay in the queue 
SC = customer - average service time 


WC = customer - average wait in the system = DC + SC 
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QM 
QO 


QL 


DV 
SV 
QV 


NC 
PS 


maximum delay experienced by any customer 
probability that a customer's delay is greater than 
Zero. 

probability that a customer's delay is greater than 


XD. SeChhe- vwolue or XDias an Anpuwt!: option.) 


maXimum queue length reached during the simulation 
probability that the queue length will be greater 
than zeror 

probability that queue length will be greater than LQ 


(the value of LQ is an input option). 


standard deviation of the delay in queue 
Standard deviation of the service times 


standard deviation of the queue Length 


total number of arrivals 


percentage not lost (losses model only) 
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V. PARAMETER INPUT 


s RUNT IME 
3 Ici 
2 
i 
1 


IDIST { 1 
KK 1 1 
RATE 38.88 38.88 14.88 
Cv 1.88 1.88 1.88 
I DSPLN | 1 


LO 1 MAXO 
xD 1.58 NPOP 
TDUAL 1 P 


ENTER 1. CONTINUE 
2. STATISTICS 


Picture 7 


Parameters can be input by the card reader or by the 
teletype terminal. Im either case, the namelist form is 
used, where the entry form is NAME = X, NAME being the 
variable title and X being the value being assigned. For 
real variables, the decimal point must be entered. Name- 


list input is terminated with a *. 


A. SAMPLE INPUT CARD 


IDEV = 2, ITYPE = 3, RATE(2) = 3.0 * 
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B. NAMELIST INPUT VIA THE GRAPHICS TERMINAL 

Lfsthis mode Of entry iS “Specitied, “one variable at 
atime as entered, followed by a carriage return. Entries 
are displayed, as entered, on the graphics screen. Termina- 
von Occurs with the entry of a * followed by a carriage 


mecurn. 


c INPUT VARIABLES 
The following conventions are employed for specifying 
the options and parameters to be used. 
MYPE = 1 gives G/G/C 
2 gives G/G/C with losses 
3 gives G/G/C with finite source 
4 gives G/G/C with feedback 
5 gives two simultaneous G/G/C queues 
RUNTIME = desired length of simulation (minutes) 
me) = number of servers for queue 1 
Het = number of servers for queue 2 (only applicable if 


TETYPE = S-is chosen) 


IDIST(I) = 1. K - ERLANG 
2. HYPEREXPONENTIAL 
3. ° DEGENERATE 
KK(I) = K in K-ERLANG (use K=1 for exponential) 
RATE(I) = Rate in number per minute 
GViCb) = Coctftictent. of variation: (specify for IDIST(I) = 2 


only. This variable is computed otherwise.) 
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macre I = 1 for‘arrivals to quewe 1 (except’ finite model) 
2 £or arrivals to queue 2 or finite population 
idle parameters 
3 for service for queue 1 (except finite model) 
4 for service for queue 2 or finite population 


service parameters 


IDSPLN (TIT) Ics ~~ PEO 

(pe BAC, 

$2 shortest Service time first 
4, Random draw 


where I = 1’ for queue I] 


2 for queue 2 


LQ Darametersan the statistic P(Q > LQ) 


XD 


parameter, in, the, statastic  P(wait.> XD) 
IDUAL = 0. Independent arrival streams 
1. Identical arrival streams 


(applicable to ITYPEs="5 “ony ) 


MAXQ = maximum queue size. (for.G/G/C-with losses only) 


-NPOP = population size (for finite source option only) 


P = completion probability (for feedback option only) 
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VI. THE PROGRAM 


The program requires the use of the XDS 9300 digital 
computer, the CI 5000 analog computer, and a graphics 
terminal using the GATED monitor. 

The analog computer wiring diagram necessary to drive 
the display sequence appears in APPENDIX A. Before 
beginning, set Delay Flop Switches DFOO, DFO1, and DFO2 
£0.12 MS and set counter FOO to 1 and FO4 to 0. Then select 
RUN and DIGITAL COMPUTER on the analog control panel. This 
Provides for movement interrupts to occur at: a rate of 
100/11 hertz. If desired, the timing values will be dis- 
played on the RATIOMETER by selecting TRUNK 420 prior to 
DIGITAL COMPUTER. 

The program consists of a FORTRAN program and a META- 
SYMBOL 9300 program. Because the history is buffered out 
to the device assigned the logical unit value 7, an addi- 
erenal control card is needed. 

Job cards ee deck run should be: 
JOB 
AAGT 
AASSIGN 7=DFI1A 
aAFORTRAN LS,GO 


FORTRAN DECK 
—~META 9300 SI,LO,GO 


META-SYMBOL DECK 


Zo 
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PROGRAM STEPS 
SET-UP ANALOG 
COMPUTER 
RESET GATED 
GRAPHICS MONITOR) 
SPECIFY 
GRAPHICS TERMINAL 
INPUT PARAMETERS INPUT PARAMETERS 
VIA TERMINAL VIA CARDS 
DISPLAY DISPLAY 
STATISTICS PARAMETERS 
DISPLAY 
STATIC EVeTURE 
SIMULATION 
CREATE HISTORY 
QUEUING SYSTEM 
DYNAMIC DISPLAY 
WRITE PARAMETERS 
AND STATISTICS 
DISPLAY DISPLAY 
PARAMETERS STATISTICS 
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APPENDIX A 


ANALOG SCHEMATIC 


SCOPE. Y¥2 
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PSET OQ &4 


T 420 
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